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Abstract 

We present a preliminary measurement of CP-violating parameters S and C from fits of the 
time-dependence of meson decays to r]'K^ . The data were recorded with the BABAR detector at 
PEP-II and correspond to 232 x 10^ BB pairs produced in e^e annihilation through the T{4S) 
resonance. By fitting the time-dependent CP asymmetry of the reconstructed B^ — > r]'K^ events, 
we find 5 = 0.60 ± 0.31 ± 0.04 and C = 0.10 ± 0.21 ± 0.03, where the first error quoted is statistical 
and the second is systematic. We also perform a combined fit using both r]'Kg and rj'K^ data, and 
find S = 0.36 ± 0.13 ± 0.03 and C = -0.16 ± 0.09 ± 0.02. 
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1 Introduction 

Decays of mesons to charmless hadronic final states such as (j)K^, K~ , r]'K^, tt^K^ and 
fo{980)K^ proceed mostly via a single penguin (loop) amplitude with the same weak phase as in 
B decays to a charmonium state plus a meson (l|,2j. However Cabibbo-Kobayashi-Maskawa 
(CKM)-suppressed amplitudes and new particles in the loop can introduce other weak phases whose 
contribution is not negligible j2 IB] • 

Fig. n^a) shows the diagram describing the B — B mixing. The amplitudes shown in Fig. E^b)- 
(d) are relevant for the decay B^ — > rj'K^. All of the amplitudes are suppressed by small CKM 
matrix elements, but the tree diagram for B^ shown in Fig. Wi^) is expected to be smaller [21 1^ 
since there is additional CKM suppression and color suppression. 




(c) (d) 

Figure 1: Feynman diagrams describing (a) B — B mixing, the decay B^ — > rj'K^ via (b, c) internal 
gluonic penguin and (d) color-suppressed tree. 

For the decay B^ rj'K^, the additional contributions of other weak phases from CKM- 
suppressed amplitudes are expected to be small, so the time-dependent CP asymmetry measurement 
provides an approximate measurement of sin2/3. We define S and C as the coefficients of sine and 
cosine oscillation terms, respectively, in the BB decay rate distributions (see Eqn. below). AS 
is the deviation between S in the decay B^ — > r]'K^ and S, equal to sin2/3, in the charmonium 

decays. Theoretical bounds for this deviation have been calculated with an SU(3) analysis 
[HI [H] . Such bounds have been improved by the measurements of B^ decays to a pair of neutral 
charmless light pseudoscalar mesons (51111) with the conclusion that AS* is expected to be less than 
0.10 (with a theoretical uncertainty less than ~0.03 due to the assumptions in the calculation). 
QCD factorization calculations conclude that AS* is even smaller ^Oj- A significantly larger value 
of AS* could arise from non-Standard- Model amplitudes 

The CP-violating asymmetry in the decay B^ — > r]'K^ has been measured previously by the 
BABAR and Belle jJH experiments using the final state tj'Kg. In the measurement presented 
in this paper we use events reconstructed in the i]'K^ final state. We present also the measurement 
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of CP-violating asymmetry combining present rj'K^ data with the rj'K'g data used in the previous 
BABML measurement jllj. 

2 The BABAR Detector and Dataset 

The results presented in this paper are based on data collected in 1999-2004 with the BABAR 
detector at the PEP-II asymmetric e^e~ collider jjl] located at the Stanford Linear Accelerator 
Center. An integrated luminosity of 211 fb^^, corresponding to about 232 million BB pairs, was 
recorded at the T{AS) resonance ( "on-resonance" , center-of-mass energy y/s = 10.58 GeV). 

The asymmetric beam configuration in the laboratory frame provides a boost of /?7 = 0.56 to 
the T{AS). Charged particles are detected and their momenta measured by the combination of 
a silicon vertex tracker (SVT), consisting of five layers of double-sided detectors, and a 40-layer 
central drift chamber, both operating in the 1.5 T magnetic field of a solenoid. The tracking system 
covers 92% of the solid angle in the center-of-mass (CM) frame. 

Charged-particle identification (PID) is provided by the average energy loss [dE/dx) in the 
tracking devices and by an internally reflecting ring- imaging Cherenkov detector (DIRC) covering 
the central region. A K/tt separation of better than four standard deviations [a) is achieved for 
momenta below 3 GeV/c, decreasing to 2.5 a at the highest momenta in the B decay final states. 
Photons and electrons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). The EMC 
provides good energy and angular resolutions for detection of photons in the range from 30 MeV to 
4 GeV. The energy and angular resolutions are 3% and 4 mrad, respectively, for a 1 GeV photon. 

The flux return (IFR) for the solenoid is composed of multiple layers of iron and resistive plate 
chambers for the identification of muons and long-lived neutral hadrons. 

3 Event Selection and Analysis Method 

Monte Carlo (MC) simulations of the signal decay modes, BB backgrounds, and detector response 
are used to establish the event selection criteria. We reconstruct r]' mesons through the decays 
rj' — > p^'^ {il'p-y) v' ~^ ??vr+7r^ {"n'-q-Kn) with r] 77. The photon energy in laboratory system 
must be greater than 50 MeV for rj candidates, and greater than 200 MeV in ij'p^. We make the 
following requirements on the invariant masses (in MeV/c^): 490 < m^^< 600 for 7], 510 < mTTTr < 
1000 for /9°, 945 < m,,^^ < 970 for rj'^^^, and 930 < nip^ < 980 for rj'^^. We require the PID 
information of the signal pions to be consistent with the pion hypothesis. 

Signal candidates are reconstructed from clusters of energy deposited in the EMC or from 
hits in the IFR not associated to any charged track in the event. Because the energy of the -fC^ 
is not measured, we determine the candidate laboratory momentum from its flight direction 
determined from the rj' vertex and the centroid of the EMC (or IFR) candidate and the constraints 
of and B^ masses to their nominal values 15_. 

A B meson candidate is characterized kinematically by the energy difference AE = E^ — ^y/s, 
where E'^ is the CM energy of the B meson. Signal events are peaked within ±10 MeV of AE = 
while background events extend towards positive values of AE (this is a consequence of the mass 
constraint used to determine the momentum). We require —0.01 < AE < 0.08 GeV. This 
choice is dictated by the need of preserving a region with enough background for a fit to that 
component. 

To reject background due to continuum e'^e" — > qq events {q = u, d, s, c), we make use of the 
angle 9t between the thrust axis of the B candidate and that of the rest of the tracks and neutral 



9 



clusters in the event, calculated in the center-of-mass frame. The distribution of cos 9x is sharply 
peaked near ±1 for combinations drawn from jet-like qq pairs and is nearly uniform for the isotropic 
B meson decays; we require | cos^rl < 0.8 in the r]' r/Tr^vr" subdecay mode and | cos 9t\ < 0.75 
in the rj' — > p^j subdecay. 

For further suppression of continuum background we require that the total missing transverse 
momentum projected along the direction, where the total momentum is calculated with all 
charged tracks and neutral clusters (without the K^, is no more than 0.45 GeV/c lower than the 
calculated transverse momentum of the K'^ candidate. We also require that the cosine of the polar 
angle of the total missing momentum in laboratory system to be less than 0.95. 

The purity of the candidates reconstructed in the EMC is further improved by a cut on 
the output of a neural network (NN) that takes cluster-shape variables as its inputs. The NN 
was trained using MC signal events and data events in the side band distribution (defined as 
0.04 < AE < 0.08 GeV). We validated the performance of the NN using candidates in the 
reconstructed — > J/ipKf^ events. 

All selection criteria have been chosen using MC signal and background events to maximize the 
expected statistical significance of signal yield in the data. 

The BB backgrounds were estimated using Monte Carlo simulations of B^B^ and B^B^ . We 
found a small evidence of BB background from b c decays in the sub-decay mode rj' — > p^j , so 
we added this component to the fit. 

For each B^ v'^l candidate (Bcp), we reconstruct the decay vertex of the other B meson 
(i?tag) from the remaining charged tracks in the event and identify its flavor. The time difference 
At = tcp — ttag, where tcp and ttag are the proper decay times of the Bqp and -Btag, respectively, 
is obtained from the measured distance between the Bqp and i?tag decay vertices and from the 
boost (/?7 = 0.56) of the e~*"e~ beam system. The distribution of At is: 

F{At) = — {1 =F Aw ± (1 - 2uj) [-rjS sm{AmdAt) - C cos(AmrfAt)]} , (1) 

4r 

where the upper (lower) sign denotes a decay accompanied by a B^ (B^) tag) t is the mean B^ 
lifetime, Am^ is the mixing frequency, ry is the CP eigenvalue of the final state (?? = -|-1 for 
rj'K^, rj = —\ for rj'K'g) and the mistag parameters uj and Aoj are the average and difference, 
respectively, of the probabilities that a true B^ {B^) meson is tagged as B^ {B^). The tagging 
algorithm, based on six tagging categories, is an improved version of what was used in the previous 
BABAR publication Separate neural networks are trained to identify primary leptons, kaons, 
soft pions from D* decays, and high-momentum charged particles from B decays. Each event is 
assigned to one of the tagging categories based on the source of tagging information and on the 
estimated mistag probability. The distribution F{At) is convolved with a resolution function to 
account for the finite vertex resolution of the detector. 

4 Maximum Likelihood Fit 

We use an unbinned, multivariate maximum-likelihood fit to extract signal yields and CP violating 
parameters. We use the following discriminating variables: AE, a Fisher discriminant J^, At. In the 
decay mode rj' we add the rj' mass and the variable 7i, defined as the cosine of the p meson's 

rest frame decay angle of a pion with respect to the tj' flight direction. The Fisher discriminant 
combines four variables: the angles with respect to the beam axis of the B momentum and B thrust 
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axis in the T(4S') frame, and the zeroth and second angular moments of the energy flow excluding 
the B candidate around the B thrust axis. 

We indicate with j the event species: signal, continuum background, or BB background. For 
each species j and each flavor-tagging category c, we define a total probability density function 
(PDF) for an events i as: 

VI, = V,{AE^) V.iD V,{Ml,)V,{W) V.iAf, a^,, c). (2) 

where is the error on At for an event i. We define the extended likelihood function for the 
input events in category c as 

( 

£c = exp I - ^ Ujfj^A n(^sig/sig>cKig,c + ^qqfqq,c'P\q + ^BBf BB,c^Bb)^ (3) 

where rij is the number of events with species j, and fj^c is the fraction of category-c events with 
species j. We fix fsig,c and f^B c fsfiav^c, the values measured with a sample of neutral B decays 
to flavor eigenstates, Bfiav 

The total likelihood function for all categories is given as the product of the likelihoods over 
the seven tagging categories (including a category for untagged events for yield determination). 



5 Results 

The reconstruction efficiency is 10.3% and 11.6% for 'q'^^K^ and ri'^^^K^ respectively. In Tabled 
we give the number of the signal yield and the parameters S and C. Note that the sign of the CP 
eigenvalue of the final state is out of the definition of S parameter (see Eqn. ©), so S parameter 
has the same sign in rj'Kg and r]'Kf^ events. The rj'K^ data are those used in BABAR previous 
measurement Combining -q'^^^K^ events and rj'p^Kl events, we measure S = 0.60 it 0.31 and 
C = 0.10 lb 0.21. In this fit we have 42 free parameters: S, C, signal yields (2), BB background 
yield (1), continuum background yields (2) and fractions (12), background At, AE, T , r]' mass and 
Ti PDF parameters (23). In the final fit, combining tj'K^ and t]'K^, we have 138 free parameters: 
S, C, signal yields (7), BB background yield (3), continuum background yields (7) and fractions 
(42), background At, AE, T , rj' mass and H. PDF parameters (77). 



Table 1: Results with statistical errors for the B^ — > r]'K^ time-dependent fits (decays with in 
upper part on the table and decays with in lower part of the table). 



Mode 


Signal yield 


S 


C 




137 ±22 
303 ± 49 


0.38 ± 0.44 
0.88 ±0.43 


0.34 ±0.29 
-0.15 ±0.29 



ri~r^+ _ 188 ±15 0.01 ±0.28 -0.18 ±0.18 

rj'K'^ . 430 ±26 0.44 ±0.19 -0.30 ±0.13 

^(3.)..^°+.- 54 ±8 0.79 ±0.47 0.11 ±0.35 

%(77K.^Xo 44 ±9 -0.04 ±0.57 -0.65 ± 0.42 

r?' K° 89 ± 23 -0.45 ± 0.68 0.41 ± 0.40 



Combined fit 1245 ± 67 0.36 ±0.13 -0.16 ±0.09 
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Figure 2: The likelihood ratio C{Sg)/[C{Sg) + C{Bg)] for ir,i,^Kl (left) and ip.,Kl (right). The 
points represent the on-resonance data, the histograms are from PDF generated events of back- 
ground (blue area) and background plus signal (shaded red area). 

The agreement between PDF simulated events and data is investigated using likelihood ratios. 
We generate signal and background Monte Carlo samples from the PDFs and, using the fitted 
parameter values from nominal fit, we calculate the likelihoods for both samples. In Fig. Elwe show 
the likelihood ratio C{Sg) /[C{Sg) + C{Bg)] for the two sub-decays rj'p^Kl and rj'^^^K^ for data and 
for the PDF generated events. In Fig. we show the projection onto /S.E of a subset of the data 
for which the signal likelihood (computed without the plotted variable) exceeds a mode-dependent 
threshold that optimizes the sensitivity. We show in Fig. |1] the At projection and asymmetry for 

6 Systematic Uncertainties and Cross-checks 

The contributions to the systematic uncertainties in S and C for rj'K^ are summarized in Table [21 

We evaluate the uncertainties associated with the PDF shapes by variation of the parameters 
describing each discriminating variable. Systematic errors associated with signal parameters (At 
resolution function, tagging fractions, and dilutions) are determined by varying their values within 
their errors. Uncertainties due to Am^^ and tb are obtained by varying these parameters by the 
uncertainty in their world average values All changes are combined in quadrature obtaining 
an error of 0.02 for 5 and 0.01 for C. 

We vary the SVT alignment parameters in the signal Monte Carlo events by the size of mis- 
alignments found in the real data, and assign the resulting shift in the fit results as the systematic 
error of 0.01 for both S and C . 

The systematic errors due to interference between the CKM-suppressed h — > ucd amplitude 
and the favored b cud for some tag-side B decays are found to be negligible for 5 and gives a 
contribution to the C uncertainty of about 0.01. 

The systematic error due to BB background is estimated to be 0.03 in S and 0.01 in C 
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Figure 3: Projection onto AE for —>■ rj' (sum of the sub-decay modes rf'^^K^ and -q'^^^K^) of 
a subset of the data for which the signal hkehhood (computed without the plotted variable) exceeds 
a mode-dependent threshold that optimizes the sensitivity. Points with errors represent the data, 
solid curve the full fit functions and dashed blue curve the total background functions. 




Figure 4: Projections onto At for r]'K^ (sum of the sub-decay modes rj'p^K^ and tj'^^^K^) of a subset 
of the data for which the signal likelihood exceeds a mode-dependent threshold that optimizes the 
sensitivity. Points with errors represent the data, solide curve the full fit functions, and dashed 
blue curve the total background functions, for (a) and (b) B^ tagged events. The asymmetry 
between B^ and B^ tags is shown in (c). 
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parameter. An uncertainty of 0.01 is assigned to account for limitations of Monte Carlo statistics 
and modeling of the signal. We assign an uncertainty of 0.01 to account for the uncertainty in the 
position and size of the beam spot, determined from variation of these quantities in signal MC. 
The total systematic error is obtained by summing individual errors in quadrature. 

Table 2: Estimates of systematic errors for r]'K^. 



Source of error 


aiS) 


a{C) 


PDF Shapes 


0.02 


0.01 


SVT alignment 


0.01 


0.01 


Tag-side interference 


0.00 


0.01 


BB Background 


0.03 


0.01 


MC statistics/modeling 


0.01 


0.01 


Beam spot 


0.01 


0.01 


Total 


0.04 


0.03 



We have also performed a number of checks of our results. When we fit the combined sample 
rj'Kg and Tj'Kf^ with the value for C fixed to zero, we find S = 0.37 it 0.13. We produce samples of 
pseudo-experiments generated with events produced to match the PDF distributions. From these 
samples, we verify that the fit bias on S and C is negligible and that there is a good agreement 
between expected and observed errors. The fit was also verified with our B^ J/t/jKf^ data sample. 

7 Conclusion 

In a sample of 232 million BB pairs we have reconstructed 137 it 22 rj'^^^Kf^ events and 303 it 49 
rj'p^K^ events. We use these events to measure the time-dependent asymmetry parameters 5* and 
C: 

S = 0.60 ± 0.31 (stat) ± 0.04 (syst) 
C = 0.10 ±0.21 (stat) ±0.03 (syst) 

Using this sample and the rj'K'g sample found in Ref. we obtain a total of 1245 ± 67 rj'K^ 
events and with a combined fit of all data we measure: 

S = 0.36 ±0.13 (stat) ±0.03 (syst) 
C = -0.16 ±0.09 (stat) ±0.02 (syst) 

All these results are preliminary. Our result for S differs from the BABAR value of sin 2/3 = 
0.722 ± 0.040 ± 0.023 in charmonium decays UH] by 2.8 standard deviation. 
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